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Theoretical and experimental study of high-acceleration patterns in vibrated granular layers is 
presented. The order parameter model based on parametric Ginzburg-Landau equation is used to 
describe strongly nonlinear excitations including hexagons, interface between flat anti-phase do- 

I — ^ ' mains and new localized objects, super-oscillons. The experiments confirmed the existence of super- 

oscillons and bound states of super-oscillons and interfaces. On the basis of order parameter model 
we predict analytically and confirm experimentally that the additional subharmonic driving results 

C^-H ' in controlled motion of the interfaces. 

o 

^^ . PACS: 47.54.+T, 47.35.-(-i,46.10.-(-z,83.70.Fn 
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The collective dynamics of granular materials is a subject of current interest R-fl]. The intrinsic dissipative nature of 
the interactions between the constituent macroscopic particles give rise to several basic properties specific to granular 
substances and setting granular matter apart from the conventional gaseous, liquid, or solid states. 

Driven granular systems manifest collective fluid-like behavior: convection, surface waves, and pattern formation 
(see e.g. [|l|). One of the most fascinating examples of this collective dynamics is the appearance of long-range 
^ ■ coherent patterns and localized excitations in vertically-vibrated thin granular layers |&fl • The particular pattern is 
'~^, determined by the interplay between driving frequency / and acceleration of the container T — An'^Af^/g {A is the 

. ' amplitude of oscillations, g is the gravity acceleration) |^,||. 
►j^ I Patterns appear at F « 2.4 almost independently of the driving frequency /. At small frequencies / < /* HQl the 

^\^ • transition is subcritical (hysteretic), leading to formation of squares. In the hysteretic region, localized excitations 
^\ such as individual oscillons and various bound states of oscillons appear as F is decreased. For higher frequencies 
^^ f > f* the onset pattern is stripes, and at the frequency slightly higher than /* the transition becomes supercritical. 

^^ Both squares and stripes, as well as oscillons, oscillate at half of the driving frequency //2. At higher acceleration 
(F > 4), stripes and squares become unstable, and hexagons appear instead. Further increase of acceleration at 
F w 4.5 converts hexagons into a domain-like structure of flat layers oscillating with frequency //2 with opposite 
phases. Depending on parameters, interfaces separating flat domains, are either smooth or "decorated" by periodic 
j^ undulations. For F > 5.7 various quarter-harmonic patters emerge. 

C^ [ The pattern formation in thin layers of granular material was studied theoretically by several groups. Direct 

I . molecular dynamics simulations ^^ (see also Q) reproduced a majority of patterns observed in experiments and 

'"O many features of the bifurcation diagram, although until now have not yielded oscillons and interfaces. Hydrodynamic 

H : and phenomenological models [0,12| reproduced certain experimental features, however neither of them offered a 

systematic description of the whole rich variety of the observed phenomena. In Refs. |13,H we introduced the order 

parameter characterizing the complex amplitude of sub-harmonic oscillations. The equations of motion following from 

the symmetry arguments and mass conservation reproduced essential phenomenology of patterns near the threshold 

of primary bifurcation: stripes, squares, and, localized objects, oscillons. 

In this paper we describe high acceleration patterns on the basis of the same order parameter model and compare it 
d ' with experimental observations. Our preliminary results were published earlier in Refs. |g,|lj,|l^. Here we show that 
at large amplitude of driving both hexagons and interfaces emerge. We find morphological instability leading to the 
formation of "decorated" interfaces. We study the motion of the interface under the influence of a small subharmonic 
component in the driving acceleration. We also find a new localized structure, " super-oscillon" , which exists for high 
acceleration values. We discuss possible mechanisms of saturation of the interface instability. Our experimental results 
demonstrate the existence of super-oscillons and bound states of super-oscillons and interfaces. They also confirm 
theoretical predictions of the externally controlled interface motion. 

The structure of the article is following. In Sec. I we introduce our phenomenological order-parameter model. In 
Sec. II we develop the weakly-nonlinear theory for the flat period-doubled state of the vibrated layer. In Sec. Ill we 
analyze the interface solution and study its stability with respect to transverse perturbations. In Sec. IV we discuss 
new localized solutions. In Sec. V we present the combined theoretical phase diagram of the model. In Sec. VI we 



■0^ caused by the spatially 

, p!4| to describe the pattern 

(which corresponds to low 



demonstrate that in a certain limit our model can be reduced to the real Swift-Hohenberg equation. This equation 
also possesses a similar variety of patterns including stripes, hexagons, stable and unstable interface solutions, as well 
as localized solutions. In Sec. VII we demonstrate that additional subharmonic driving results in drift of the interface 
with the velocity determined by the amplitude of the driving and the direction determined by the relative phase. 
In Sec. VIII we present our experimental results. It includes the study phase diagram and the effect of additional 
subharmonic driving. Section IX summarizes our conclusions. 

I. PARAMETRIC GINZBURG-LANDAU EQUATION 

The essence of the model |13| , |l4| is the order parameter equation for the complex amplitude '0 of parametric layer 
oscillations h — ip exp(i7r/i) + c.c. at the frequency //2 coupled with the equation for the thickness of the layer p 
averaged over the period of vibrations, 

dt^jj = 70* - (1 - iLu)^P + (1 + i5)VV - IV'I V - P0 (1) 

dtp = CV-{pV\ij\')+pV^p (2) 

Here 7 is the normalized amplitude of forcing at the resonant frequency /. Linear terms in Eq. (n^) can be obtained 
from the complex growth rate for infinitesimal periodic layer perturbations h ~ exp[A{k)t + ikx]. Expanding A(fc) for 
small k, and keeping only two leading terms in the expansion A(fc) — — Aq — Aifc^ we obtain the linear terms in Eq. (Q), 
where b — ImAi/ReAi characterizes ratio of dispersion to diffusion and parameter to = (Oq — 7r/)/i?eAo, Qq — —ImAo 
characterizes the frequency of the driving. In Eq. (0) C, /? are transport coefficients for p. Slowly- varying thickness of 
the layer p controls the dissipation rate (the last term in Eq.(|l])). The second equation (|) describes re-distribution 
of the averaged thickness due to the diffusive flux ex — Vp, and an additional flux ex — pV 
nonuniform vibrations of the granular material. This coupled model was used in Refs. [|l3,| 
selection near the threshold of the primary bifurcation. It was shown that at small CPo/3 
frequencies and thick layers), the primary bifurcation is subcritical and leads to the emergence of square patterns. 
For higher frequencies and/or thinner layers, transition is supercritical and leads to roll patterns. At intermediate 
frequencies (/ « /*), the stable localized solutions of Eqs.(|l|),(||) corresponding to isolated oscillons and variety of 
bound states were found in agreement with experiment. 

In this paper we focus on high-acceleration patterns at high frequencies. In Ref. p3[ it was indicated that the 
density transport coefficient /3 is proportional to the energy of plate vibration (oc /^, and, therefore, it should 
increase with the driving frequency /. As a result, for high frequencies the coupling between p and tp becomes less 
relevant, and one can assume p =const and exclude it from Eq.(n^) by rescaling. Then the model can be reduced to a 
single order-parameter equation 

dt^ = 7-0* - (1 - itj)0 + (1 + i6)VV - 1012^; (3) 

Eq. (^ also describes the evolution of the order parameter for the parametric instability in vertically oscillating fluid 
layers (see 0,0). 

It is convenient to shift the phase of the complex order parameter via ip — 0exp(i0) with sin 20 ~ lo/^. The 
equations for real and imaginary part ^j — A + iB are: 

dtA = {s - I) A - 2ojB - {A^ + B^)A + \/^{A- hB) (4) 

dtB = -{s + l)B - {A^ + B^)B + V\B + hA) (5) 



II. STABILITY OF UNIFORM STATES 

At s < 1, Eqs. (Q),(||) has only one trivial uniform state A — 0, i? = 0, At s > 1, two new uniform states appear, 
A — ±^/s — 1, B = 0. The onset of these states corresponds to the period doubling of the layer flights sequence, 
observed in experiments ^ and predicted by the siniplc inelastic ball model |^,|l^ • Signs ± reflect two relative phases 
of layer flights with respect to container vibrations ||l^ . 

The trivial flat state A = B = loses stability if the following condition fulfilled (compare fl^ ) : 



Because of the symmetry of Eq.(||), small perturbations from the trivial state lead to the formation of a periodic 
sequence of rolls or stripes. 

Let us analyze the stability of the non-trivial states A — ±\/s — 1, B = with respect to small perturbations with 
wavenumber k 

At the threshold the uniform state loses its stability with respect to periodic modulations with the critical wavenumber 
kc at s < Sc (correspondingly, 7 < 7c), where 



'^ = 26 ^^) 

1 + 52 ^^' 

Small perturbations with every direction of the wavevector grow at the same rate. The resultant selected pattern 
is determined by the nonlinear competition between the modes. In the presence of the reflection symmetry ^ -^ —ip, 
quadratic nonlinearity is absent, and cubic nonlinearity near the trivial state favors stripes corresponding to a single 
mode. Near the fixed points A = zLy/s — 1,B = the reflection symmetry for perturbations U ^ —U, y^— V^is 
broken, and hexagons emerge at the threshold of the instability. To clarify this point we perform weakly-nonlinear 
analysis of Eqs. (^.(||) for s = Sc — e, and e ^ 1. At e ^ 0, the variables U and V are related as in linear system: 

* = ^ Aj exp[ikr] -I- c.c, (10) 

where |A;| = kc and ry = [2{sc - 1) + fc^] /(6fc2 - 2uj). 
Corresponding adjoint eigenvector is 



L/n /I 
V+J \r]+ 



(11) 



where 77+ = — [2(sc — 1) + k^] /bk^- Substituting Eq. ( |lO|) into Eqs. (|4|),(^ and performing the orthogonalization, 
we obtain equations for the slowly- varying complex amplitudes Aj, j = 1,2,3 (we assume only three waves with 
triangular symmetry, favored by quadratic nonlinearity) 

dtAj = 2eAj + a2A*+iA*_i 

-a3(|A,f + 2(|A,_i|2-H|A,+i|2))A„ (12) 

where the coefficients 02, 03 are 



I + rj- 



02 = ±2V7^T 2 + ^ , ' ,03 = 3(1 + 7777+). (13) 

V l + riri+ J 

Eqs. dll) are well-studied (see [pOJ). There are three critical values of e: ea = —02/4003, en = 02/203, and 
eB = 2a^Ja3. The hexagons are stable for e^ < e < es, and the stripes are stable for e > eu- Thus, near s = Sc 
the model exhibit stable hexagons ||21[] . Since we have two symmetric fixed points, both up- and down- hexagons 
co-exist. For smaller s stripes are stable, and for larger s, flat layers are stable, in agreement with observations P,f5|, 
as well as direct numerical simulations of Eq.(y) ||22|] . The above analysis requires the values of eA,B.R to be small. 
For parameters uj,b — 0{\), this requirement is satisfled for e^, but not for €b.r- The estimates can be improved 
by substituting s = Sc + e instead of s = Sc in ( [l3| ) . The resulting range of stable hexagons is plotted in the phase 
diagram (^,7) (see below Fig.^. 

III. INTERFACE SOLUTION 

At s > 1, Eqs. (Q),(|3) have an interface solution connecting two uniform states A = ii/s — 1, B = 0. Let x be the 
coordinate perpendicular to the interface and y is parallel to the interface. The interface is the stationary solution to 
Eqs. i),!!): 



1)Aq - 2ujBo - (Al + B^)Ao + Aq^-^ - bBo^^ = 0, 
-{s + l)Bo - {Al + Bl)Bo + Bo^a, + bAo^a, = 0. 



(14) 



For 6 = Eqs. (|lj) have a solution of the form ^o = ±a/s — ltanh(-\/s — la;/2), Bq = 0. For b ^ the solution is 
available only numerically. We used shooting matching technique to find this solution in the range of parameters w, b 
and 7. A typical solution is shown in Fig. H^. As one sees from the Figure, for b ^ the asymptotic behavior of the 
interface exhibits decaying oscillations. 
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FIG. 1. Interface solution to Eqs. y_^ obtained numerically for to — l,b — 4,^ — 2. Solid line corresponds to A and dashed 
line to B, respectively. 



Let us now consider the perturbed solution 



Ao 
Bo 



Aix) 
B{x) 



exp[A(A:)t + iky]. 



(15) 



where k is the transverse modulation wavenumber and \{k) is the corresponding growth rate. For A^ B we obtain 
linear equations 



L 



B 



{X{k) + k^) 



bk^ 



B 

~A 



(16) 



where the matrix L is of the form 

L = 



3^2 -Bl + dl, -2uj - 2AqBq - bdl 
-2A0B0 + bdl --5 - 1 - ^0 - 3Sg + dl 



(17) 



In order to determine the spectrum of eigenvalues X{k) we have solved Eq. (|l^) along with stationary Eqs. (|lj) 
numerically using numerical matching-shooting technique. We have found that for small uj the interface is stable with 
respect to transverse undulations, i.e. A(fc) < 0. However, for the values of lu above certain critical value tod'y, b) the 
interface exhibits transverse instability: A(fc) > in the band of wavenumbers |fc| < kc, see Fig. 0. 
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FIG. 2. Spectrum of eigenvalues A(A;) for lo = 1.3, 6 = 4 and 7 = 2 (solid line) and a; = 1, 6 = 4 and 7 = 2 (dashed line) 
obtained from numerical solution of Eqs. [16l 

This instability is confirmed by direct numerical simulations of Eq.(|3|). An example of the evolution of slightly 
perturbed interface is shown in Fig. H. Small perturbations grow to form a "decorated" interface. With time these 
decorations evolve slowly and eventually form a labyrinthine pattern. 

The neutral curve for this instability can be determined as follows. Numerical analysis shows that at the threshold 
the most unstable wavenumber is A: = and we can expect that for A: — > A ~ fc^, see Fig. g. Expanding Eqs.([l6|) in 
power series of fc^ : 






Ad) 

5(1) 



(18) 



in the zeroth order in k we obtain L(A(°), B'") = 0. The corresponding solution is the translation mode A^^^ 
dxAfi{x\ B^^> = dxBo{x). In the first order in k^ we arrive at the linear inhomogeneous problem 



S(i) 



(A(fc) 



^ \b(0) 



tk^i ^"^ 



-A(o) 



(19) 



A bounded solution to Eq. (jlj) exists if the r.h.s. is orthogonal to the localized mode of the adjoint operator A+, B^ 
The adjoint operator L+ is of the form 



L^ 



l-3Al-Bi+ dl -2AoBo + bdl 
-2uj - 2AoBo - bdl, -s-l-Al-iBl + dl 



(20) 



Since the operator L"*" is not self-adjoint, the adjoint mode does not coincide with the translation mode, and, therefore, 
must be obtained numerically. 

The orthogonality condition fixes the relation between A and k: 



A = -crfc^,cr = 1 + 6 



J-^(A(0)5+_S(0)A+)^^ 

/_!°^(A(o)A+ + B(o)B+)da; 



(21) 



where a is "surface tension" of the interface. The instability, corresponding to the negative surface tension of the 
interface, onsets at cr = 0. Although the interface solution A{x) is not localized (see Fig. |l|), the integrals in Eq. ( |2l| ) 
converge because the adjoint functions A^,B^ arc localized. The numerically-obtained adjoint eigen-mode is shown 
in Fig. 1^. The neutral curve is shown in Fig. |j. Direct numerical simulations of Eq.(g) show that this instability 
leads to so called "decorated" interfaces (see Fig. H,a), however At later stages the undulations grow and finally form 
labyrinthine patterns (Fig. |j, b-d). 




FIG. 3. Adjoint eigenfunction obtained numerically for uj = 1.3, fo = 4, 7 = 2. Solid line corresponds to A^ and dashed line 
to B^ , respectively. 



FIG. 4. (a-d) Interface instability and labyrinth formation, Eq. 
snapshot are taken at times t = 1000, 1600, 2300, 4640. 




ctJ = 2, 6 = 4, 7 = 2.9, domain size 100 x 100 units, the 



The emergence of the labyrinthine patterns as a result of interface instabihty contradicts experiments |Jj5[| where 
stable decorated interfaces are typically observed. Thus, the model Eq. (0) does not capture saturation of the interface 
instability and requires modification. We also checked that the coupling to the density field introduced in Ref. [l^ ] 
also docs not provide desired saturation. 

Let us now discuss possible mechanisms of saturation of the transversal instability of the interface, which is not 
captured by Eq. (H). The reason for proliferation of the labyrinths is that local stabilization mechanisms cannot 
saturate the "negative surface tension instability" of the interface. Indeed, the behavior of small perturbation C, of 
almost flat interface is described by the following linear equation: Q — — crAC, where a is the surface tension. The 
linear growth must be counter-balanced by nonlinear terms. Due to translation invariance the local nonlinear! ty can 
be function of gradients of ^ only. In the lowest order the first term is |VCp, since the quadratic term does not 
saturate monotonic instability. Since for the modulation wavenumber fc ^ one has |VCP ^ fc^lCP and A^ ~ fc^C, 
the linear instability always overcomes local nonlinearity at small wavenumbers. As a result small perturbations grow 
and saturation happens due to non-local interactions as in a labyrinthine pattern which is stabilized by non-local 
repulsion of the fronts (see Q ) . 

As we see, in our model the interface instability does not result in stable decorated interfaces. One of the possibility 
for the saturation of this instability could be local convection induced by the interface. Indeed, due to shear motion 
of grains at the interface (flat regions moves in anti-phase) local convective flow can be created. The scale of this flow 
can be large then the interface thickness. One of the possible manifestation of the convection beside grain transport 
is reduction of the effective "viscosity" of fiat layer, or increase in the sensitivity to external vibrations [p3. As we 



can see from Fig. [^ increase in vibration magnitude 7 suppresses the interface instability. As the instabihty develops, 
the length of the interface grows, which increases the convection and finally suppresses the instability. This feedback 
mechanism results in saturation of the interface instability and creates stable decorations. 

Although the details of this coupling are not known yet, it can be implemented within the order-parameter approach 
in a phcnomcnological way. In Eq. the forcing 7'0* must be replaced by (70 + w)tp* , where w includes the effect of 
local convective flow and 70 = const. We take the simplest possible form of the coupling to the convective flow: 



/ dr' exp 



.'|2 



l^Hr'W 



(22) 



where rg characterizes the typical scale of the convective flow and e is the amplitude of the coupling, and S is the 
area of integration domain. Since |VV'(?'')I i^ nonzero only at the interface, this integral is proportional to the total 
length of the interface. With the growth of the interface length the value of "effective forcing" 7 = 70 + w grows as 
well, and eventually leaves the instability domain. 

We performed numerical simulations with Eq. (g) modified in this fashion. The convolution integral (^) was 
calculated using the fast Fourier transform. The results are presented in Fig . H. Remarkably, for certain range 
of parameters the instability indeed saturates due to nonlocal term Eq. (122), and we find interfaces with stable 
decorations, see Fig. ^,a. 




FIG. 5. Computer simulations of modified Eq. (H). Surface plots show Re^ for different stationary solutions: (a) the 
saturated interface, uu = 1.2,7 = 1-85, b — 4, e — 0.002, ro — 40, L = 160, (b) saturated interface with bound super- oscillons, 
uj = 1.2,7 = 1.75,6 = 4, e = 0.002, ro = 40, L = 160. (c) single super- oscillon aX uj ^ 1.2,7 = 2.2, fo = 4, e = 0.002, ro = 40, 
L = 100 . 

IV. LOCALIZED SOLUTIONS 

In addition to the interface solution which exists for s > 1, Eq. ^ for even higher values of s possesses a 
variety of other nontrivial two-dimensional solutions including localized " super-oscillons" . The typical solution and 
its corresponding localized linear mode are shown in Fig. g. 
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FIG. 6. (a) Axi-symmetric localized solution (super-oscillon) obtained for b = 4,uj = 1 and 7 = 1.8 for Eq. 
^-solid line, B- dashed line, (b) Localized eigenfunction corresponding to the maximum eigenvalue A = —0.0009. 



at 6 = 4, 



Although super-oscillon appears similar to the "conventional oscillon" discovered in Ref. Q, there is a significant 
difference. Asymptotic behavior of the super-oscillon is ^ — > ^0 7^ for r ^ 00 whereas the regular oscillon has 
^ ^ 0. As a result super-oscillons always have the same phase as their background. Oppositely-phased super- 
oscillons must be separated by the interface connecting two opposite phases of the period-doubled flat layer whereas 
oppositely-phased oscillons can form stable bound states (dipoles, H). 

Super-oscillons exist in a relatively narrow domain of parameters of Eq. (|^) . Decrease of 7 leads to destabilization 
of the localized eigen-mode leading to subsequent nucleation of new " super-oscillons" on the periphery of the unstable 
super-oscillon, with the increase of 7 super-oscillons disappear in a saddle-node bifurcation. 

Interestingly, for certain values of 70 we also find an interesting structure consisting of a decorated interface and 



a chain of super-oscillons bound to the decorations, see Fig. |^,b. These objects can exist completely independently 
of the interface, Fig. g c, as one may expect for Eq. (g). Thus, we conclude that extended Eq. (y) qualitatively 
describes observed phenomenology. 



V. PHASE DIAGRAM 



The results of linear stability analysis can be summarized in the phase diagram shown in Fig. 0. Remarkably, 
the structure of the phase diagram F igJ^ is qualitatively similar to that of the experiments for high frequencies of 
vibration, see Fig. |^ and Refs. p-|5|,p|,|15f. Increasing vibration amplitude 7 leads to transition from a trivial state to 
stripes, hexagons, decorated interfaces, and finally, to stable interfaces. Transition from unstable to stable interfaces 
also occurs with decreasing to (increasing vibration frequency /), in agreement with Ref. [ [l5[|2^ . In experiments, at 
yet higher F, quarter-harmonic patterns appear, however these patterns are not described by our model. Our analysis 
also predicts co-existence of stripes and hexagons in a very narrow parameter range (see Sec. II). Super-oscillons were 
found in a narrow domain close to line 7 of Fig. [7|. 

Note that square patterns which are emerge at lower frequencies, are not included in the phase diagram. They are 
obtained within the full model Eqs.(|l|),(||) with additional density field p, see fl^ . 



^-2 




FIG. 7. Phase diagram for Eq. (|3|) at 6 = 4. Line 1, 7^ = (ij + 6)^/(1 + b^)\ line 2, a^ = Y - ^ = 1; line 3, s = s^; line 4, 
a = Sc — £«; line 5, s = Sc — es; line 6, s = Sc — ea; line 7, surface tension a = G. Below line 1 there are no patterns, between 
lines 1 and 4 stripes are stable, between lines 5 and 6 hexagons are stable, above line 2, non-trivial flat states exist, and above 
line 3 they are stable. Interfaces are unstable below line 7. Hexagons and stripes coexist between lines 4 and 5 

VI. SWIFT-HOHENBERG EQUATION 

Near the line s — 1 (see Fig. R) Eq. (H) can be simplified. In the vicinity of this line A ^ [s — 1)^/^ and 
B r^ [s- if/"^ < A. In the leading order we can obtain from Eq. (||) B = 6V^A/2, and Eq. (|) yields [|§: 



dtA ={s- I) A -A^ + {1- ujh)V'^A V^A. 



(23) 



Rescaling the variables t — > (s — l)t, A —^ (s — 1) ^/^A, x — > (2(s — l)/6-^)^/*, we arrive to the Swift-Hohenberg 
equation (SHE) 



dtA = A-A^-5\/'A- VM 



(24) 



where 



Lob- 1 



s-1 



(25) 



The description by SHE is valid ii S ^ 0(1) which implies additional condition for the validity of this approach 
w6- 1 < 1 at s ^ 1. 

Although this equation is simpler than Eq. (0), it captures many essential features of the original system dynamics, 
including existence and stability of stripes and hexagons in different parameter regions (see P7| , p8|) , existence of the 
interface solutions, interface instability, super-oscillons and emergence of labyrinthine patterns ||29|. Indeed, simple 
analysis shows that the growth rate of the instability of the uniform state A = 1 as a function of the perturbation 
wavenumber is determined by the formula X{k) = —2 + 6k^ — k'^, so it becomes unstable at S > Sc = 2-\/2 at critical 
wavenumber kc = \/2- As in the original model, near the threshold of this instability, subcritical hexagonal patterns 
are preferred. Interface stability can also be analyzed more simply as the linearized operator corresponding to the 
model (24) is self-adjoint. The threshold value of 5 is obtained from the following solvability condition 



(S 



th^Ox ~ '^^Oxx) dx — 



(26) 



which yields S^^ — 1.011. It yields the equation for the neutral curve of the interface instability in the Swift- Hohcnberg 
limit in the form: 



1 



2{ujb- 1) 



(27) 



Fig. g, a-d shows the development of the interface instability within SHE with subsequent transition to labyrinthine 
patterns, similar to the dynamics of the full parametric equation (0). Again, to saturate the instability, an additional 
non-local mechanism is required. 
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FIG. 8. (a-d) labyrinth formation from a circular spot, SHE, 5 = 1.4, domain size 100 x 100 units, t = 200, 1300, 1700, 1900. 



As it was shown recently, Swift-Hohenberg equation also possesses localized solutions (see [^9|j30|). These localized 
states are analogous to the super-oscillons of Eq. (g) . 



VII. SUBHARMONIC FORCING AND INTERFACE MOTION 

Let us focus on the effect of additional subharmonic driving on the motion of the interface. In the original model 
Eq. (^ the interface does not move due to symmetry between left and right halves of the interface. However, if the 
plate oscillates with two frequencies, / and //2, the symmetry between these two states connected by the interface, 
is broken, and interface moves. The velocity of interface motion depends on the relative phase of the sub-harmonic 
forcing with respect to the forcing at /. The effect of a small external subharmonic driving applied on the back-ground 
of primary harmonic driving with the frequency / can be described by the following equation (cf. Eq. (ph) 



^tV = 7^-* - (1 - Jtj)V' + (1 + i6)VV ■ 



■0 + ge' 



(28) 



where q characterizes the amplitude of the sub-harmonic pumping, and $ determines its relative phase. 

For small g, we look for moving interface solution in the form tp = ipt^ix — vt) -\- qipi{x — vt) -I- ... and v = 0{q), 
alternatively. 
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Bq{x — vi) 



a{x — vt) 
b[x — vt) 



Oiq') 



(29) 



Solvability condition (compare with Eq. (|2l|)) fixes the interface velocity as a function of amplitude and the phase 
of external subharmonic driving 



cos $ / A+da; + sin <f> / B+dx 
' J{A+d,Ao + B+d,Bo)dx 



gasin($ ~ $0) 



(30) 



where a = const has the meaning of susceptibility of the interface for external forcing and $0 — const is some phase 
shift. The explicit answer is possible to obtain for & = when A^ = dxAo, and $ = 0,7r, which yields the interface 
velocity v = ^^qAq^ = ^^q{s—l)~^. In general, A, B, ^+, i3+, and hence v, can be found numerically. The interface 
velocity as function of q, $ is shown in Fig. |9[ 
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FIG. 9. Interface velocity v for lj = 1, 6 = 4, 7 = 2.5 vs g at $ = 0. Inset: w vs $ at g = 0.01. (•) - numerical results, ( — ) - 
analytical expression (pd). 

Thus, from this analysis we conclude that additional subharmonic driving results in controlled motion of the 
interface. The velocity of the interface depends on the amplitude of the driving, and the direction is determined by 
the relative phase $. 

VIII. EXPERIMENTAL RESULTS 

We performed experiments with a thin layer of granular material subjected to periodic driving. Our experimental 
setup is similar to that of Q. We used 0.15mm diameter bronze or copper balls, and the layer thickness in our 
experiments was 10 particles. We performed experiments in a rectangular cell of 4 x 12 cm^ . We varied the acceleration 
F and the frequency / of the primary driving signal as well as acceleration 7, frequency /i and the phase $ of the 
secondary (additional) driving signal. The interface position and vertical accelerations were acquired using a high- 
speed video camera and accelerometers and further analyzed on a Pentium computer. To maintain and measure the 
proper acceleration at / and //2 wc employ the lock-in technique with our signal originating from an accelerometer 
attached to the bottom of the cell. This allows for the simultaneous real-time feed back and control of the device. To 
reduce the interstitial gas effects |3l|| we reduced the pressure to 2 mTorr. Our visual data is acquired using a high 
speed digital camera (Kodak SR-lOOOc), in addition to high speed recording this also allows for a synchronization 
between the patterns and the image acquisition. 

A. Phase Diagram 

The experimental phase diagram is shown in Fig. |l^. It is similar to that of Ref. |^,^, however the transition 
from hexagons to interfaces is elaborated in more detail. The dashed line indicates the stability line for interfaces 
with respect to periodic undulations. To determine this stability limit, we kept the value of acceleration F fixed and 
increased the driving frequency /. For each value of the frequency we extracted the interface width by averaging up 
to 10 images of the interface. In order to find that amplitude of the periodic undulations I we subtracted from the 



11 



obtained value the thickness of the flat interface Iq. The resultant amplitude or undulations / as function of driving 
frequency / is shown in Fig. |ll|. As one sees from the Figure, the amplitude of undulations / decreases gradually 
approaching the critical value of the frequency. However, a small hysteresis at the transition point cannot be excluded 
because of large error bars. 
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FIG. 10. Phase diagram. Material is 150 nm copper balls, layer thickness 10 particle diameters. 
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FIG. 11. The amplitude of periodic undulations Z as a function of driving frequency /. The acceleration amplitude F — 4.1 



As it follows from the phase diagram, for small amplitude of the vertical acceleration F < 3.6 stripes are the only 
stable pattern (note that we focus on high-frequency patterns when squares do not exist, compare |3|). At slightly 
higher F = 3. 6.. 3. 7 stripes and hexagons coexist. For higher F — 3. 7.. 4.0 the hexagons becomes stable. Due to the 
sub- harmonic character of motion both up and down hexagons may co-exist separated by a line phase defect, Fig. 
n3,a. There exists a narrow band, from F = 3.85 to F = 4.0, where localized states, or super-oscillons, appear both 
within the bulk of the material and also pinned to the the front between the interface and the bulk Fig. |lj,b. For 
even higher values of the acceleration F > 4.0 the super-oscillons disappear (Fig. |l^,c) leading to isolated interfaces 
with periodic decorations. 
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FIG. 12. Representative high-acceleration patterns in the rectangular A x 12 cm cell, driving frequency / = 40 Hz: (a) 
r = 3.75 shows the onset of the interface with up and down hexagons, (b) V = 3.94 super- oscillons are present in the flat layer 
and are pinned to the front near the interface, (c) F = 4, isolated decorated interface. 

B. Controlled Motion of Interface 

In the absence of additional driving the interface drifts toward the middle of the cell (see Fig. O). We attribute 
this effect to the feedback between the oscillating granular layer and the plate vibrations due to the finite ratio of 
the mass of the granular material to the mass of the vibrating plate. Even in the absence of subharmonic drive, the 
vibrating cell can acquire subharmonic motion from the periodic impacts of the granular layer on the bottom plate 
at half the driving frequency. If the interface is located in the middle of the cell, the masses of material on both 
sides of the interface are equal, and due to the anti-phase character of the layer motion on both sides, an additional 
subharmonic driving is not produced. The displacement of the interface X from the center of the cell leads to a mass 
difference Am on opposite sides of the interface which in turn causes an additional subharmonic driving proportional 
to Atti. In a rectangular cell Atti cx X. Our experiments show that the interface moves in such a way to decrease the 
subharmonic response, and the feedback provides an additional term ~X/t in the r.h.s. of Eq. (pO), yielding 



dX 
~dt 



'X/t + qa sin(<f> - $0) 



(31) 



The relaxation time constant r depends on the mass ratio (this also holds for the circular cell if X is small compared 
to the cell radius). Thus, in the absence of an additional subharmonic drive {q — 0), the interface will eventually 
divide the cell into two regions of equal area (see Fig. [l^ ). 

In order to verify the prediction of Eq. (^ij) we performed the following experiment. We positioned the interface 
off center by applying an additional subharmonic drive. Then we turned off the subharmonic drive and immediately 
measured the amplitude fi of the plate acceleration at the subharmonic frequency [Q . The results are presented in 
Fig. O. The subharmonic acceleration of the cell decreases exponentially as the interface propagates to the center 
of the cell. The measured relaxation time r of the subharmonic acceleration increases with the mass ratio of the 
granular layer and of the cell with all other parameters fixed. The mass of the granular layer was varied by using two 
different cell sizes: circular, diameter 15.3 cm, and rectangular, 4x12 cm, while keeping the thickness of the layer 
unchanged. For these cells we found that the relaxation time r in the rectangular cell is about four times greater than 
for the circular cell (see Fig. |l^). This is consistent with the ratio of the total masses of granular material (52 grams 
in rectangular cell and 198 grams in circular one). In a separate experiment the mass of the cell was changed by 
attaching an additional weight of 250 grams to the moving shaft, which weights 2300 grams. This led to an increase 
of the corresponding relaxation times of 15-25 %. The relaxation time r increases rapidly with F (see Fig. [l3| inset). 



13 



0.02 



zL 0.01 



0.00 




5 10 

t (sec) 



FIG. 13. Amplitude ol the subharmonic acceleration /i of the cell averaged over 4 measurements vs time for the interface 
propagating to the center of the cell for F — 3.97 and / = 40 Hz. Circles/squares correspond to the circular/rectangular cells, 
open/closed symbols correspond to light/heavy cells, respectively. Heavy cells differ from light ones by an additional weight of 
250 g attached to the moving shaft. Solid lines show exponential fit /x ^ exp(i/r) + const. Inset: r vs F for light rectangular 
ceU. 



When an additional subharmonic driving is applied, the interface is displaced from the middle of the cell. For small 
amplitude of the subharmonic driving 7, the stationary interface position is X = qar sin($ — $0)7 since the restoring 
force balances the external driving force. The equilibrium position X as function of <f> is shown in Fig. |l4|a. The 
solid line depicts the sinusoidal fit predicted by the theory. Because of the finite mass ratio effect, the amplitude of 
the measured plate acceleration ^ at frequency /i also shows periodic behavior with $, (see Fig. 14 d), enabling us 
to infer the interface displacement from the acceleration measurements. For even larger amplitude of subharmonic 
driving (more than 4-5 % of the primary driving) extended patterns (hexagons) re-emerge throughout the cell. 

The velocity Vq = qa which the interface would have in an infinite system, can be found from the measurements of 
the relaxation time r and maximum displacement X^ at a given amplitude of subharmonic acceleration 7, Vq = X,„/r, 
see Eq.(^l|). We verified that in the rectangular cell the displacement X depends linearly on 7 almost up to values at 
which the interface disappears at the short side wall of the cell (see Fig. |lj, inset). Figure |l^ shows the susceptibility 
a as a function of the amplitude of the primary acceleration F. The susceptibility decreases with F. The cusp-like 
features in the F-dependence of a (and t, see inset to Fig. 03), are presumably related to the commensurability 
between the lateral size of the cell and the wavelength of the interface decorations. 
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FIG. 14. (a) Equilibrium position X and (b) amplitude of measured subharmonic acceleration fi as functions of phase $. 
Circular cell, T = 4.1, / = 40 Hz, g = 1.25% of T. 
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FIG. 15. Susceptibility a vs F at / = 40Hz, rectangular cell. Inset: Displacement X as function of g at $ = 260". 
We developed an alternative experimental technique which allowed us to measure simultaneously the relaxation 
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time r and the "asymptotic" velocity Vq. This was achieved by a small detuning A/ of the additional frequency /i 
from the exact subharmonic frequency //2, i.e. A/ = /i — //2 <C /. It is equivalent to the linear increase of phase 
shift $ with the rate 27rA/. Th is linear growth of the phase results in a periodic motion of the interface with frequency 
A/ and amplitude X,„ = VqJ \Jt~'^ + (27rA/)2 (see Eq (pl|)). The measurements of the "response functions" X„i(A/) 
are presented in Fig. n6l From the dependence of Xm on A/ we can extract parameters Vq, ol and t by a fit to the 
theoretical function. The measurements are in a very good agreement with previous independent measurements of 
relaxation time r and susceptibility a. For comparison with the previous results, we indicate the values for r and a, 
obtained from the response function measurements of Figs. O and hS (stars). The measurements agree within 5 %. 
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FIG. 16. Maximum displacement Xm from center of rectangular cell as function of frequency difference A/ = /i — //2 for 
/ = 40 Hz, and for V = 3.97 (circles) and V = 4.1 (diamondsV Dashed lines are fit to X™ = VqI^Jt-"^ + (27rA/)2. The values 
of a and r obtained from the fit are also indicated in Figs. FL3 and [13 by arrows with stars. 



IX. CONCLUSIONS 



In this paper we studied the dynamics of thin vibrated granular layers in the framework of phcnomenological model 
based on the Ginzburg-Landau type equation for the order parameter characterizing the amplitude of sub-harmonic 
sand vibrations. Although rigorous derivation of this model is not possible up to date, our approach gives a significant 
insight in the problem in hand and results in testable predictions. Our model is complementary to more elaborate fluid 
dynamics-like continuum descriptions or direct molecular dynamics simulations and allows us to carry out analytical 
calculations of various regimes and predict their stability. A challenging problem is to establish a quantitative relation 
between our order parameter model and other models and experiments. 

We have shown in this paper that the generic parametric Ginzburg-Landau model Eq. (^ captures not only patterns 
of vibrated sand near the primary bifurcation, but also large acceleration patterns such as hexagons, interfaces, 
and super-oscillons. The structure of the phase diagram FigJ^ for high frequencies and amplitudes of vibrations 
is qualitatively similar to that of previous experiments Refs. P,pp|p^, as well as the experiments reported here. 
Increasing vibration amplitude leads to transition from a trivial state to stripes, hexagons, decorated interfaces, and 
finally, to stable interfaces. Interface decorations were found to be a result of the transversal instability of a flat 
interface which is analogous to the negative surface tension. We found that Eq. (|h is not sufficient to describe 
the saturation of interface instability. We propose a possible non-local mechanism of saturation of this instability 
which takes into account the dependence of the overall length of the interface and magnitude of parametric forcing. 
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which may occur due to a large scale convective flow induced by the sand motion near the interface. We described 
analytically the motion of the interface under the symmetry-breaking influence of the small subharmonic driving. 

In our experimental study we found that on a qualitative level the theoretical phase diagram is similar to the 
experimental one. The experiments also confirmed existence of super-oscillons and their bound states with the interface 
as it was predicted by the theory. Further experimental study is necessary to elucidate the specific mechanism for 
saturation of the interface instability. 

Stimulated by the theoretical prediction, we also performed experimental studies of interface motion under addi- 
tional subharmonic driving. The experiment confirmed that the direction and magnitude of the interface displacement 
depend sensitively on the amplitude and the relative phase of the subharmonic driving. Moreover, we found that 
period-doubling motion of the flat layers produces subharmonic driving because of the finite ratio of the mass of the 
granular layer and the cell. This in turn leads to the restoring force driving the interface towards the middle of the 
cell. 
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